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Abstract. Existence of 2D enstrophy cascade in a suitable mathematical setting, 
and under suitable conditions compatible with 2D turbulence phenomenology, is 
^-^ known both in the Fourier and in the physical scales. The goal of this paper is to 

show that the same geometric condition preventing the formation of singularities - 
i -Holder coherence of the vorticity direction - coupled with a suitable condition 
on a modified Kraichnan scale, and under a certain modulation assumption on 
f\] evolution of the vorticity, leads to existence of 3D enstrophy cascade in physical 

scales of the flow. 
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,-H ! l. Introduction 

The vorticity-velocity formulation of the 3D Navier-Stokes equations (NSE) reads 

(1.1) d t u + (u ■ V)w = (w • V)w + Aw 

where u is the velocity and u = curl u is the vorticity of the fluid (the viscosity is 
set to 1). A comprehensive introduction to the mathematical study of the vorticity 
can be found in 1MB02I . 

I/") ■ In 2D, the vortex-stretching term, (u> ■ V)u, is identically zero, and the nonlin- 

earity is simply a part of the transport of the vorticity by the flow. 

r * This allowed the authors to adopt to 2D a general mathematical setting for the 

f—^ ■ study of turbulent cascades and locality in physical scales of 3D incompressible 

viscous and inviscid flows introduced in |DaGrll-l| and |D aGrll^2| , respectively, 
to establish existence of the enstrophy cascade and locality in 2D [Da Gr 11-31 . 

Since the vortex-stretching term is not a flux-type term, the only way to estab- 
lish existence of the enstrophy cascade in 3D in this setting is to show that its con- 
5_j ■ tribution to the ensemble averaging process can be suitably interpolated between 

integral-scale averages of the enstrophy and the enstrophy dissipation rate. 

This is where coherent vortex structures come into play A role of the coher- 
ent vortex structures in turbulent flows was recognized as early as the 1500's in 
Leonardo da Vinci's "deluge" drawings. On the other hand, Kolmogorov's K41 
phenomenology IKo41-l[ IKo41-2 1 does not discern geometric structures; the K41 
eddies are essentially amorphous. As stated by Frisch in his book Turbulence, The 
Legacy of A.N. Kolmogorov |Fr95|, "Half a century after Kolmogorov's work on 
the statistical theory of fully developed turbulence, we still wonder how his work 
can be reconciled with Leonardo's half a millennium old drawings of eddy motion 
in the study for the elimination of rapids in the river Arno." This remark was fol- 
lowed by a discussion on dynamical role, as well as statistical signature of vortex 
filaments in turbulent flows. Several (by now classical) directions in the study of 
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the vortex dynamics of turbulence - both in 2D and 3D - are presented in Chorin's 
book Vorticity and Turbulence (cf. [Ch94| and the references therein). The approach 
exposed in |Ch94] is essentially discrete (probabilistic lattice models); on the other 
hand, the first rigorous continuous statistical theory of vortex filaments was given 
by P.-L. Lions and Majda in ILMOOI . 

Local anisotropic behavior of the enstrophy i.e., self-organization of the re- 
gions of high vorticity in coherent vortex structures - most notably vortex fila- 
ments/tubes - is ubiquitous. In particular, local alignment or anti-alignment of the 
vorticity direction, i.e. local coherence, is prominently featured in turbulent flows. 
A strong numerical evidence, as well as several theoretical arguments explaining 
the physical mechanism behind the formation of coherent structures - including 
rigorous estimates on the flow directly from the 3D NSE - can be found, e.g., in 
ll^9^ [SJO9j lCPS95[lGGH97llGFD99llOh09llGMlll . One way to look at the phe- 
nomenon of local coherence of the vorticity direction is to interpret it as a mani- 
festation of the general observation that the regions of high fluid intensity are - in 
the vorticity formulation - locally 'quasi 2D-like' (in 2D, the vorticity direction is 
globally parallel or antiparallel). 

The rigorous study of geometric depletion of the nonlinearity in the 3D NSE was pi- 
oneered by Constantin when he derived a singular integral representation for the 
stretching factor in the evolution of the vorticity magnitude featuring a geomet- 
ric kernel depleted by coherence of the vorticity direction (cf. BCo9 4|). This was 
followed by the paper [CoFe93| where Constantin and Fefferman showed that 
as long as the vorticity direction in the regions of intense vorticity is Lipschitz- 
coherent, no finite-time blow up can occur, and later by the paper [daVeigaBe02| 
where Beirao da Veiga and Berselli scaled the coherence strength needed to deplete 
the nonlinearity down to \ -Holder. 

A full spatiotemporal localization of the vorticity formulation of the 3D NSE was 
developed by one of the authors and the collaborators in [ GrZh06 . Gr09 . GrGulO^Tl 
GrGulO-2J. The main obstacle to the full localization of the evolution of the enstro- 
phy was the spatial localization of the vortex-stretching term, (cj • V)u. An explicit 
local representation for the vortex-stretching term was given in [Gr09J, and the 
leading order term reads 

_&_ 1 

iBixo^Rj^^dxidyk \x- 

where e/fcj is the Levi-Civita symbol and cf) is a spatiotemporal cut-off associated 
with the ball B(x a ,R). The key feature of dl.2b is that it displays both analytic (via 
a local non-homogeneous Div-Curl Lemma |GrGulO-2[) and geometric (via coher- 
ence of the vorticity direction [Gr09 GrG ulO-ll ) cancelations, inducing analytic and 
geometric local depletion of the nonlinearity in the 3D vorticity model. (For a different 
approach to localization of the vorticity-velocity formulation see [ChKaLe07].) 

The present work is envisioned as a contribution to the effort of understanding 
the role that the geometry of the flow and in particular, coherent vortex structures, 
plays in the theory of turbulent cascades. More precisely, we show - utilizing the 
aforementioned localization within the general mathematical framework for the 
study of turbulent cascades in physical scales of incompressible flows introduced 
in [ Da Gr 11-11 (in this case, via suitable ensemble-averaging of the local enstrophy 
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equality) - that coherence of the vorticity direction, coupled with a suitable condition 
on a modified Kraichnan scale, and under a certain modulation assumption on 
evolution of the vorticity, leads to existence of 3D enstrophy cascade in physical 
scales of the flow. This furnishes a mathematical evidence that, in contrast to 3D 
energy cascade, 3D enstrophy cascade is locally anisotropic, providing a form of a 
reconciliation between Leonardo's and Kolmogorov's views on turbulence on the 
enstrophy level. 

It is worth pointing out that our theory of turbulent cascades in physical scales 
of 3D incompressible flows is on the energy level |DaGrll-l, DaGrll-2, Da Grl2-l| 
consistent both with the K41 theory of turbulence and the Onsager's predictions 
on existence of the inviscid energy cascade (and consequently, with the phenomena 
of anomalous dissipation and dissipation anomaly), as well as with the previous rig- 
orous mathematical work on existence of the energy cascade in the wavenumbers 
[FMRT01J. In particular, on the energy level, it does not 'see' geometric structures. 
It is only here, i.e., on the enstrophy level, that the role of coherent vortex struc- 
tures is revealed. A distinctive feature of our theory that makes incorporating the 
geometry of the flow and in particular, local coherence possible is the fact that the 
cascade takes place in the actual physical scales of the flow. To the best of our 
knowledge, there is nothing in the K41 theory that would contradict existence of 
the 3D enstrophy cascade; on the other hand, given that K41 takes place primar- 
ily in the Fourier space, i.e., in the wavenumbers, formulating conditions (within 
K41) faithfully reflecting various geometric properties of the flow is a much more 
challenging enterprise. 

The paper is organized as follows. Section 2 recalls the ensemble-averaging 
process introduced in [DaGrll-lJ, and Section 3 the spatiotemporal localization of 
the evolution of the enstrophy developed in |GrZh06IIGr09llGrGulO-lUGrGulO-2ll . 
Existence and locality of anisotropic 3D enstrophy cascade is presented in Section 

4. 



2. Ensemble averages 

In studying a PDE model, a natural way of actualizing a concept of scale is to 
measure distributional derivatives of a quantity with respect to the scale. Let xq 
be in B(0, Rq) (Rq being the integral scale, B(0, 2R ) contained in ft where Q is the 
global spatial domain) and < R < Rq. Considering a locally integrable physical 
density of interest / on a ball of radius 2R, B(xq, 2R), a local physical scale R - 
associated to the point xo - is realized via bounds on distributional derivatives 
of / where a test function ip is a refined - smooth, non-negative, equal to 1 on 
B(xq, R) and featuring optimal bounds on the derivatives over the outer i?-layer 
- cut-off function on B(xq, 2R). More explicitly, 

(2.1) \(D a f,i;)\ < f l/HD-VI < (c(«)il/U 5(Q) ) 

for some c(a) > and 5(a) in (0, 1). (This is reminiscent of Bernstein inequalities 
in the Littlewood-Paley decomposition of a tempered distribution.) 
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Henceforth, we utilize refined spatiotemporal cut-off functions <p = 4> x „.r.t = 
ipj], where n = rfr{t) G C°°(0,T) and ip = ip Xo , R (x) e T>(B(xq,2R)) satisfying 

(2.2) < 77 < 1, ?7 = Oon(0,T/3), ry = 1 on (2T/3,T), ^<^ 
and 

(2.3) 0<*<1, * = 1« B (,, B ), M<|, J^<§, 

for some | < pi,/»2 < 1. 

In particular, <f> = ipof) where i/jq is the spatial cut-off (as above) corresponding 
to s = and R = R . 

For x near the boundary of the integral domain, S(0,Rq), we assume addi- 
tional conditions, 

(2.4) < i> < il> 

and, if B{x ,R) <£ B(0,R ), then V € V{B{0 7 2R )) with ip = 1 on B{x ,R) n 
B(0, Ro) satisfying, in addition to l!2.3l l, the following: 

ip = i/jq on the part of the cone centered at zero and passing through 

5(0, Rq) n B(x , R) between 5(0, R ) and 5(0, 2Rq) 



(2.5) 
and 



ip = on 5(0, i?o) \ B(xq 7 2R) and outside the part of the cone 
(2.6) centered at zero and passing through 5(0, Ro) n B(x ,2R) 

between 5(0, Rq) and 5(0, 2R ). 

A physical scale R - associated to the integral domain B(0, Rq) - is realized via 
suitable ensemble-averaging of the localized quantities with respect to \K\, K 2 )~ 
covers' at scale R. 

Let K\ and K 2 be two positive integers, and < R < Rq- A cover {B(xi, R)}™ =1 
of the integral domain B(0, R ) is a (K\ , K 2 )-cover at scale R if 

-Rj ^ n£Kl {-R 

and any point x in B(Q, R ) is covered by at most K 2 balls B(xi,2R). The param- 
eters K\ and K 2 represent the maximal global and local multiplicities, respectively. 

For a physical density of interest /, consider time-averaged, per unit mass - 
spatially localized to the cover elements B(xi, R) - local quantities j Xi ,R, 

fx it R= 7ri -ml f(x,t)<p s RT (x,t)dxdt 

1 JO K JB(x t ,2R) 

for some < 6 < 1, and denote by (F)r the ensemble average given by 

1 " A 

(f)r = - y Jxj,r- 

i=l 

The key feature of the ensemble averages {(F)r}o<r<r is that (F)r being sta- 
ble, i.e., nearly-independent on a particular choice of the cover (with the fixed pa- 
rameters K\ and K 2 ), indicates there are no significant fluctuations of the sign 
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of the density / at scales comparable or greater than R. On the other hand, if / 
does exhibit significant sign-fluctuations at scales comparable or greater than R, 
suitable rearrangements of the cover elements up to the maximal multiplicities - 
emphasizing first the positive and then the negative parts of the function - will 
result in (F) r experiencing a wide range of values, from positive through zero to 
negative, respectively. 

Consequently, for an a priori sign-varying density, the ensemble averaging pro- 
cess acts as a coarse detector of the sign-fluctuations at scale R. (The larger the maximal 
multiplicities K\ and K 2 , the finer detection.) 

As expected, for a non-negative density /, all the averages are comparable to 
each other throughout the full range of scales R, < R < Ro; in particular, they 
are all comparable to the simple average over the integral domain. More precisely, 



(2.7) 


-^F < (F) R < K*F 
A* 


for all < R < R , where 




F = 


f J ^3 J f(x,t)4(x,t 



andK* = K*{K 1 ,K 2 ) > 1. 

There are several properties of the ensemble averaging process that although 
being plausible, deserve a precise analytic description/quantification. Besides the 
above statement on non-negative densities, perhaps the most elemental one is that 
if the averages at a certain scale are nearly independent of a particular choice of 
a (Ki, K 2 ) -cover (a"i and K 2 fixed), then essentially the same universality property 
should propagate to larger scales. In order to obtain a precise quantitative prop- 
agation result, the universality is assumed on an initial interval, rather than at a 
single scale (this is due to the presence of smooth cut-offs with prescribed rates of 
change). Two types of results are presently available. 

Let / be a locally integrable function (a density), and K\ and K 2 two positive 
integers. 

TYPE I. Assume that there exists i?» > such that for any R in [R*, 2i?„) and 
any (27 Ki, l^A^-cover at scale R, the averages (F)r are all comparable to some 

value F„; more precisely, — -F* < {F)r < C\F*. Then, for all R > 2R* and all 

(Ki, K 2 )-covers at scale R, the averages (F)r satisfy — -F* < (F)r < C 2 F*. 

TYPE II. Assume that there exists i?* > such that for any R in {\R* , ^R*) and 

1 



any (K u K 2 )-cover at scale R, -^-F* < (F) R < CiF». Then, for all R > §i?* and 
all (K\, K 2 ) -covers at scale R, the averages (F) R satisfy 



5 (f) F,<(F )R <C 3 Q F,. 

Shortly, in a Type I result, the universality is assumed with respect to more re- 
fined covers, while in a Type II result, the non-exactness of the propagation caused 
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Figure 1. 



f(x) = cos 2 (a; + 5) sinU(x - l) 2 ) 



by the smooth cut-offs is reflected in a correction to the universal value F* by the 
ratio of the scales R and R* . 

The proofs are quite long and technical, and will be provided in a separate 
publication together with computational results describing the general statistics 
of the variation of the ensemble averages of multi-scale sign-fluctuating densi- 
ties across the scales. Here, we provide a sample computation of the ensemble 
averages of a time-independent (to emphasize the spatial behavior) ID density 

f(x) = cos 2 (a: + 5) sin(i(a: - l) 2 ) (Figure 1), for R = 10 and K x = K 2 = 3; its 

global average is approximately —0.003880. The y-values in Figure 2 represent the 
ensemble averages with respect to the {K\, ii^-covers exhibiting maximal pos- 
itive and negative bias, across the range of scales from 10 -2 to 10 1 = i?o (the 
x-values represent the powers of 10); the red line slightly below the a>axis corre- 
sponds to the value of the global average. The response of the ensemble averages 
to sign-fluctuations at several different scales is clearly visible. 



3. Spatiotemporal localization of evolution of the enstrophy 

The localization of the vorticity formulation introduced in [GrZh06, Gr09| was 
performed on an arbitrarily small parabolic cylinder below a point (xo,io) con- 
tained in the spatiotemporal domain Vl x (0, T), with an eye on formulating the 
conditions preventing singularity formation at (xo,4o)- Here - in order to coor- 
dinate the notation with Section 2 - the localization will be performed on a spa- 
tiotemporal cylinder B(xq, 2R) x (0, T) where xo belongs to the integral domain 
£(0,i? )and0 <R< R . 

Suppose that the solution is smooth in B(xq, 2R) x (0, T). Multiplying the equa- 
tions by 4> to (</> = tprj being the cut-off introduced in Section 2) and integrating over 
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FIGURE 2. Ensemble averages with positive and negative bias. 



B(x ,2R) x (0, t) for some 2T/3 <t<T yields 



-|w(ir,t)| 2 V(aOda: + 



(3.1) 



IVwl 4> dx ds 



-\ui\ 2 ((j)t + A0) dx ds 



/ -|w| 2 (u • V(/>) dx ds + J (w ■ V)w • 0w dx ds. 



Suppressing the time variable, the localized vortex-stretching term can be writ- 
ten as (cf. HGr§0) 

r) 

(w • V)w • 4>uj (x) — (J)? (x) — — Uj (x) </> 5 (x) uJi (x) ujj (x) 

OXi 



-cP.V. 



d 2 1 
e jki o — « — i 1 4>^ uidy <j)* (x) UJi (x) Uj (x) + LOT 



(x ,2r) 9a;i(9y fe \x - y\ 



= -cP.V. / (u(x)xu(y))'G u (x,y)<f>i{y)^{x)dy+ LOT 

•/-B ( x ,2r) 

(3.2) = VST + LOT 

where ejki is the Levi-Civita symbol, 



(G u (x,y)) t 



d 2 



1 



k dxidy k \x - y\ 
and LOT denotes the lower order terms. 



Ui(x) 



The above representation formula for the leading order vortex-stretching term 
VST features both analytic and geometric cancelations. 

The geometric cancelations were utilized in |Gr09] to obtain a full localization 
of \ -Holder coherence of the vorticity direction regularity criterion, and then in 
BGrGulO-ll to introduce a family of scaling-invariant regularity classes featuring a 
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precise balance between coherence of the vorticity direction and spatiotemporal 
integrability of the vorticity magnitude. Denote by £ the vorticity direction, and 
let (x, t) be a spatiotemporal point, r > and < 7 < 1. A 7-Holder measure of 
coherence of the vorticity direction at (x, t) is then given by 

, + , \sm<p(i(x,t),Z(y,t))\ 

p lt r(X,t) = SUp 



I l-y 

y£B(x,r),y^x \ x V\ 

The following regularity class - a scaling-invariant improvement of | -Holder co- 
herence - is included, 

(3.3) / / \oj(x,t)\ 2 pi 2R (x,t)dxdt < 00. 

Jt -(2R) 2 JB(x ,2R) 2 ' 

On the other hand, the analytic cancelations were utilized in |Gr GulO-2| via a 
local non-homogeneous Div-Curl Lemma to obtain a full spatiotemporal localiza- 
tion of Kozono-Taniuchi generalization of Beale-Kato-Majda regularity criterion; 
namely, the time-integrability of the BMO norm of the vorticity. 



4. 3D ENSTROPHY CASCADE 

Let 1Z be a region contained in the global spatial domain f2. The inward enstro- 
phy flux through the boundary of the region is given by 

-\uj\ 2 (u ■ n) da = — / (u ■ V)oj • oj dx 
an 2 Jk 

where n denotes the outward normal (taking into account incompressibility of the 
flow). Localization of the evolution of the enstrophy to cylinder B(xq, 2R) x (0, T) 
(cf . Section 3) leads to the following version of the enstrophy flux, 

(4.1) / -\u\ 2 (u ■ V (j)) dx = - (u-V)u -(f>ujdx 

(again, taking into account the incompressibility; here, 4> = 4> Xo ,r,t defined in 
Section 2). Since V0 = (\/ip)rj, and ip can be constructed such that V-0 points 
inward - toward xq - (|4.1|l represents local inward enstrophy flux, at scale R (more 
precisely, through the layer S(x 0: R, 2R)) around the point x . In the case the point 
.To is close to the boundary of the integral domain B(0,Rq), V-0 is not exactly 
radial, but still points inward. 

Consider a (Ki,K 2 ) cover {B(xi, -R)}" =1 at scale R, for some < R < i?o- Local 
inward enstrophy fluxes, at scale R, associated to the cover elements B(xi, R), are 
then given by 

(4.2) [±\oj\ 2 (u-V<f>i)dx, 

for 1 < i < n (cj)i = <fi Xi .R.T)- Assuming smoothness, the identity Il3.lt written for 
B(xi, R) yields the following expression for time-integrated local fluxes, 
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— \u!\ 2 (u-\7(j)i)dxds= / — |w(x, t)| tpi(x) dx + j j |Vw| 2 <^ dx ds 



-|w| 2 ((0 4 ) s + A0») dxds 



(4.3) — I (ui • V)w • 4>i to dx ds 



o 
for any t in (2T/3, T) and 1 < i < n. 

Denoting the time-averaged local fluxes per unit mass associated to the cover 
element B(x i: R) by $ Xi ,R, 

(4.4) K,R = \J j^/^MVVMds, 

the main quantity of interest is the ensemble average of {&xi,R}? = i in the sense of 
Section 2; namely 

1 " 
(4-5) ($)fl = -V$, l ,i ? . 

?i * — ' 

i=l 

Since the flux density — (u • V)w • w, is an a priori sign-varying density the sta- 
bility i.e., the near-constancy of ($)# - while the ensemble averages are being run 
over all (if 1,^2) covers at scale R - will indicate there are no significant sign- 
fluctuations of the flux density at the scales comparable or greater than R. 

The main goal of this section is to formulate a set of physically reasonable con- 
ditions on the flow in -6(0, 2Rq) x (0, T) implying the positivity and near-constancy 
of ( < &)ij across a suitable range of scales - existence of the enstrophy cascade. 

We will consider the case il = M 3 . The main reason is that in the case of a 
domain with the boundary it is necessary to utilize the full spatial localization 
(joint in x and y) of the vortex-stretching term given by (|3.2|l ; this introduces a 
number of the lower order terms which - in turn - lead to terms of the form 

Cl Ry // \ UJ \ 2 $ P ~ 1(ixds 

for some 7 > 2 (c 7 is a suitable dimensional constant - scaling like i? 7 ~ 2 ). Such 
terms introduce a correction to the dissipation cut-off in the enstrophy cascade; 

namely the modified Kraichnan scale do (see (A2) below) would have to be re- 

2. 

placed by a^ (times a dimensional constant). On the other hand, in IR 3 , the full 
(spatial) localization of the vortex-stretching term can be replaced by the spatial 
localization in x only, 

(4.6) (w • V)u • 0iW (a;) = -cP.V.J(u(x) x oj(y)) ■ G^{x,y) <f>f(x) ^(x) dy; 

the y integral is then split in suitable small and large scales (similarly to [GrZh06J) 
without introducing correction terms (cf . the proof of the main result in this sec- 
tion). 

(Al) Coherence Assumption 
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Denote the vorticity direction field by £, and let M > (large). Assume that there 
exists a positive constant C\ such that 

\sm^(x,t)^{y,t))\ < Ci\x - y\^ 

for any (x, y, t) in (B(0,2R ) x B(Q,2Rq + R$) x (0,T)) n {|Vu| > M] (ip{z 1 ,z 2 ) 
denotes the angle between the vectors z\ and z-x). Shortly, I -Holder coherence in 
the region of intense fluid activity (large gradients). 

Note that the previous local regularity results HGrZh061 lGr091 imply that - un- 
der (Al) - the a priori weak solution in view is in fact smooth inside B(0, 2Rq) x 
(0, T), and can, moreover, be smoothly continued (locally-in-space) past t = T; in 
particular, we can write H.31 with t = T. 



Let us briefly remark that in the aforementioned works on regularity, the region 
of intense fluid activity is usually defined as {|w| > M} rather that {|Vu| > M}. 
Cutting-off at |w| = M here would eventually lead to replacing Eq in the definition 
of the modified Kraichnan scale (To (see the next paragraph) by 

Eo = ^J-^Jl\Vu\ 2 (b 2 Q p - 1 dxdt, 

and we prefer keeping ctq solely in terms of ui. 

(A2) Modified Kraichnan Scale 

Denote by E time-averaged enstrophy per unit mass associated with the integral 
domain 5(0, 2R ) x (0,T), 

I /' 1 /' I . ,2 ,2p-l 



Eo = Tj RJJ 2 M ^ dxdt ' 
by P a modified time-averaged palinstrophy per unit mass, 



(the modification is due to the shape of the temporal cut-off -q), and by <jq a corre- 
sponding modified Kraichnan scale, 

Then, the assumption (A2) is simply a requirement that the modified Kraichnan 
scale associated with the integral domain B(0, 2Rq) x (0, T) be dominated by the 
integral scale, 

(T < PRo, 

for a constant /?, < (3 < 1, f3 = (3(p, Ki,K 2 ,M, B T ), where B T — sup ||w(t)|| L i; 

t e(o,T) 

this is finite provided the initial vorticity is a finite Radon measure [Co90J. 
(A3) Localization and Modulation 
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The general set up considered is one of the weak Leray solutions satisfying (Al). 
As already mentioned, (Al) implies smoothness; however, the control on regularity- 
type norms is only local. On the other hand, the energy inequality on the global 
spatiotemporal domain R 3 x (0, T) implies 

tT , 

/ \uj\ 2 dxdt < oo; 
/o Jr 3 

consequently, for a given constant C 2 > 0, there exists Rq > (Rq < mm{VT, 1}; 
this is mainly for convenience) such that 

(4.7) f [ 2 \u\ 2 dxdt< i- 

for any < i?o < -^o- This is the localization assumption on i?o; the precise value 
of the constant C^ is given in the proof of the theorem - right after the inequality 

632). 

The modulation assumption on the evolution of local enstrophy on (0, T) - con- 
sistent with the choice of the temporal cut-off 77 - reads 

|u;(a;,T)| tpo(x)dx>- sup / |u>(a;,£)| i 4>o(x)dx. 
2 te(o,T) J 

Remark 4.1. Assumption (Al) is simply a quantification of the degree of local 
coherence of the vorticity direction - a manifestation of the local 'quasi 2D-like' 
behavior of turbulent flows - needed to sufficiently deplete the nonlinearity. As- 
sumption (A2) is a slight modification of the condition implying existence of the 
enstrophy cascade in 2D (cf. [DaGrll-2[); namely, a requirement that the modified 
Kraichnan scale be dominated by the integral scale. Once unraveled, it postulates 
that the region of interest exhibits large vorticity gradients (relative to the vorticity 
magnitude, in the spatiotemporal average) - indicating high spatial complexity 
of the flow - and is the enstrophy analogue of the condition implying existence 
of the energy cascade in 3D (|DaGrll-l|), i.e., large velocity gradients (relative to 
the velocity magnitude, in the spatiotemporal average). The last assumption is 
somewhat technical; however, its purpose within our theory is physical - to pre- 
vent uncontrolled temporal fluctuations of i?o-scale enstrophy. Such fluctuations 
would inevitably prevent the cascade over B(0, i?o) (the region of interest). 

Theorem 4.1. Let ubea Leray solution on R 3 x (0, T) with the initial vorticity cj being 
a finite Radon measure. Suppose that u satisfies (A1)-(A3) on the spatiotemporal integral 

domain B(0, 2R Q + i?f) x (0, T). Then, 

-^Po < (*>r < 4K* P 
for all R, jj<7 < R < R (K* > 1 is the constant in (12. 7P ). 

Proof. Throughout the proof, {K\, K2) cover parameters p, K\ and K2, as well as 
the gradient cut-off M will be fixed. Henceforth, the quantities depending only 
on p, K\ , Ki , M, Bt will be considered constants and denoted by a generic K (that 
may change from line to line). 
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Let < R < Ro- As already noted, we can write the expression for a time- 
integrated local flux corresponding to the cover element B(xi,R), JA.3J , for t = T, 



— \u)\ 2 (u ■ \7(f>i) dx ds — j —\u>(x,T)\ 2 ipi(x) dx + / / \Voj\ 2 (f>i dx ds 

t r x 

-M 2 ((0i) s + A&) dxds 
o J z 

(4.8) - / / (w • V)u ■ (f>iUJ dx ds. 



The last two terms on the right-hand side need to be estimated. 

For the first term, the properties of the cut-off <j>i together with the condition 
T>Rl>R 2 yield 

(4.9) j j±\u\ 2 {{<t>i), + &4>i)dxd8<K^j J\io\ 2 <fir l dxds. 



For the second term, the vortex-stretching term 

i-T r 

(w • V)m • (pi uj dx ds, 

the integration is first split into the regions in which | Vw| < M and | Vw| > M. In 
the first region, the integral is simply dominated by 

(4.10) K-LJ J M 2 ^- 1 dx ds 

(R < Rq < 1). In the second region, and for a fixed x, we divide the domain of 
integration in the representation formula (4.61 1, 

(lj ■ V)u • fcu (x) =-cP.V. j (uj(x) xw(y))-G u (x, y) # (x) <f>? (x) dy, 

into the regions outside and inside the sphere {y : \x — y\ = Rs}. In the first case, 
the integral is bounded by 

\u\ dy 4>i\(jj\ 2 dxds 



{|Vm|>M} J{y: \x-y\>R,i} F — V\ 

<K±mp\\u(t)\\i>[ J 'mV?'" 1 dxds 





(4.ii) - K h\ J\"\ 2 <t> 2 i p ~ ldxds - 
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In the second case - utilizing (Al), the Hardy-Littlewood-Sobolev and the Gagliardo- 
Nirenberg interpolation inequalities, and the localization part of (A3) - the follow- 
ing string of bounds transpires. 



1 ■ ■ 



\u\ dy (j>i\u)\ dxds 

{\Vu\>M} J{y: \x-y\<Rt} \X-y\2 

<K f |M| , i MfwfWsds 

JO L 2 {B{x ,2R +R§)) 2 

<K f |M| , i J|# W || 2 ||V(#u;)|| 2 d S 

JO L 2 (^B(a:o,2i? +ii d )J 

<*(/ IMP ( l\ ds Y(l SIl P IIV-J«III+/ HV(#u;)||!<fe 

\J0 L 2 (B(x ,2_Ro+fl Q 3 )J / V 2 te(0,T) Jo 

1 ^ " * ,2 , /■ T „^,I..M,2 



<g^5 U te T T W w lla + y l|V(0»||^ S 



(4.12) 



where /l"* is the constant in j2.7|l (in the last line, we used the localization assump- 
tion l|4.7t with C2 = QAK 2 K^). Incorporating the estimate 

\V(4>fu)\ 2 dx 

,2^ , „ (/N.lV u2 



< 2 / {Vuj^^dx + c / f i_pi j | w | 2 cte 

< 2 / |Vw| 2 ^ dx + K-^ j \lo\ 2 ^ 1 dx 
(4.13) 

in l|4.12|l . we arrive at the final bound, 

M dy <f>i\uj\ dxds 



{\Vu\>M}J{y.\x-y\<R3} \X - y\ 

-4^(^ te T T \\4"\\l + [^mld^+KJpj^^dx. 

(4.14) 

Collecting the estimates J4.9H4.14ll and using the modulation part of (A3), the 
relation ( 14.81 1 yields 

-\u>\ 2 (u ■ \/ (f>i) dx ds — -\uj(x,T)\ 2 Tpi(x) dx + I \\7uj\ 2 <pi dx ds + Z 

o J 2 J 2 Jo J 

(4.15) 
where 

Z < ^ (l\\^(M-,T)\\j + £ \\Vu;4\\lds^ + K^ J H 2 ^- 1 dx. 



14 R. DASCALIUC AND Z. GRUJIC 

Taking the ensemble averages and exploiting l|2.7|l multiple times, we arrive at 

1 -Po < (®)r < 4,K*P Q 



4K* 

for all icr < R < Rq, and a suitable P = 0(p,Ki,K 2 , M, B T ). D 

Remark 4.2. The first mathematical result on existence of 2D enstrophy cascade 
is in the paper by Foias, Jolly Manley and Rosa [FJMR02|; the general setting is 
the one of infinite-time averages in the space-periodic case, and the cascade is in 
the Fourier space, i.e., in the wavenumbers. A recent work [DaGrll-3| provides 
existence of 2D enstrophy cascade in the physical space utilizing the general math- 
ematical setting for the study of turbulent cascade in physical scales introduced in 
BDaGrll-ll . To the best of our knowledge, the present paper is the first rigorous 
result concerning existence of the enstrophy cascade in 3D. 

The second theorem concerns locality of the flux. According to turbulence phe- 
nomenology, the average flux at scale R - throughout the inertial range - is sup- 
posed to be well-correlated only with the average fluxes at nearby scales. In par- 
ticular, the locality along the dyadic scale is expected to propagate exponentially. 

Denoting the time-averaged local fluxes associated to the cover element B(xi,R) 
by$ Xi ,R, 

(4.16) K,R = ^J j\M\u-V<t>i)dx, 
the (time and ensemble) averaged flux is given by 

n 

(4.17) (*) fl = -V£ XiiH = J R 3 <$) R . 



The following locality result is a simple consequence of the universality of the 
cascade of the time and ensemble-averaged local fluxes per unit mass ($)r obtained 
in Theorem 4.1. 

Theorem 4.2. Let ubea Leray solution on R 3 x (0, T) with the initial vorticity uj a being 
a finite Radon measure. Suppose that u satisfies (A1)-(A3) on the spatiotemporal integral 

2 

domain B(0, 2Rq + Rq) x (0, T), and let R and r be two scales within the inertial range 
delineated in Theorem 4.1. Then 

1 /r\ 3 (tf) r 2 {r^ 3 



16K*\RJ ~ (*)* " *\R, 
In particular, ifr = 2 k Rfor some integer k, i.e., through the dyadic scale, 

1 _« . (* 



2 3fc < \*/2*R < 16R 2 2 3fe_ 



16^ 2 " (*) 

Remark 4.3. Previous locality results include locality of the flux via a smooth fil- 
tering approach presented in |E05] (see also |EA09|), and locality of the flux in 
the Littlewood-Paley setting obtained in [CCFS08J. The aforementioned results 
are derived independently of existence of the inertial range, and are essentially 
kinematic upper bounds on the localized flux in terms of a suitable physical quantity 
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localized to the nearby scales; the corresponding lower bounds hold assuming sat- 
uration of certain inequalities consistent with the turbulent behavior. In contrast, 
our result is derived dynamically as a direct consequence of existence of the turbu- 
lent cascade in view, and features comparable upper and lower bounds throughout the 
inertial range. 
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